Negative differential resistivity in superconductors with periodic arrays of pinning 

sites 
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We study theoretically the effects of heating on the magnetic flux moving in superconductors 
with a periodic array of pinning sites (PAPS). The voltage-current characteristic (V7-curve) of 
superconductors with a PAPS includes a region with negative differential resistivity (NDR) of S- 
type (i.e., S-shaped F/-curve), while the heating of the superconductor by moving flux lines produces 
NDR of N-type (i.e. , with an N-shaped V7-curve) . We analyze the instability of the uniform flux flow 
corresponding to different parts of the V7-curve with NDR. Especially, we focus on the appearance 
of the filamentary instability that corresponds to an S-type NDR, which is extremely unusual for 
superconductors. We argue that the simultaneous existence of NDR of both N- and S-type gives rise 
to the appearance of self-organized two-dimensional dynamical structures in the flux flow mode. We 
study the effect of the pinning site positional disorder on the NDR and show that moderate disorder 
does not change the predicted results, while strong disorder completely suppresses the S-type NDR. 

PACS numbers: 74.25.Qt 



I. INTRODUCTION 



Negative Differential Resistivity (NDR) and Conduc- 
tivity (NDC) can be observed in various non-linear me- 
dia. To illustrate the counterintuitive nature of this 
phenomenon, let us consider a force acting on a set of 
moving particles: NDC corresponds to a lower velocity 
of motion for these particles when the force applied to 
them increases. Two different types of NDR can be ob- 
served in the voltage-current characteristics ( V7-curves) 
of non-linear media [H, i, i, i, H, H, ■ NDR of S-type 
is characterized by the existence of three different val- 
ues of the current / corresponding to a single value of 
the voltage V. The corresponding V7-curve is S-shaped. 
A V7-curve with three different values of the voltage 
for a single value of the current is referred to as NDR 
of N-type. The corresponding VT-curve is N-shaped. 
NDR (or NDC) is commonly observed in semiconduc- 
tors, plasmas, superconductors and is used in many non- 
linear devices (see, e.g. Refs. H, i, i, i, i, S H). In 
particular, semiconductors with NDR are the basic ele- 
ments of Gunn-effect diodes and pnpn-j unctions [l|, |3]. 
V7-curves with NDR can only be observed under spe- 
cific conditions. For example, to study N-type NDR one 
has to include the corresponding sample in an electric 
circuit with fixed voltage V. Vice versa, to observe S- 
type NDR the sample should be included in a circuit 
with fixed current /. If these conditions are not ful- 
filled, the uniform current flow becomes unstable, and 
non-uniform self-organized structures (e.g., filaments and 
overheated domains with higher or lower electric fields) 
arise in the sample. Such structures are commonly ob- 
served in plasmas, semiconductors, and superconductors 



(see, e.g. Refs. Ill, 
from Refs. [H, [2,1 



5, a0)- Table 1 (using results 
1 i, la [13, [H [H, [H compares 



NDR in different non-linear media. Note, the macro- 
scopic manifestations of NDR in different media could 
be rather similar while the intrinsic physical mechanisms 
giving rise to NDR could be different. 

As a non-linear medium with NDR, we study a su- 
perconductor with artificial pinning sites. The magnetic 
flux behavior in such superconductors has attracted con- 
siderable attention due to the possibility of constructing 
samples with enhanced pinning as well as with novel and 
unusual voltage-current characteristics [I4lll5lll6lll7l.ll8l. 
113] . Present-day technology allows the fabrication super- 
conductors with well-deflned periodic arrays of pinning 
sites (PAPS). Such structures include many thousands 
of elements with controlled microscopic pinning param- 
eters. Increased interest on these systems has arisen in 
recent years, and a number of intriguing features related 
to PAPS has been revealed. 

In Ref. [20 the existence of several dynamical vortex 
phases was predicted for square PAPS subjected to per- 
pendicular magnetic field, B, close to the first matching 
field Bd, = ^o/a^, where $0 is the magnetic flux quantum 
and a is the PAPS period. The geometry of the problem 
is shown in Fig. 1 which is discussed in Sec. II. The effect 
of the dynamic phases on the VT-curve is illustrated in 
Fig. 2a. Let us assume that the field B is slightly higher 
than B^ and the number of the vortices in the sample Ny 
is higher than the number of the pinning sites Np. Let 
us now slowly increase the applied current j in the sam- 
ple. For very low current density j (phase I in Fig. 2) 
all vortices are pinned and their average velocity v is 
zero. With increasing the current j, interstitial vortices. 
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TABLE I: Comparison between non-uniform non-equilibrium states in superconductors, semiconductors, plasmas, and mangan- 
ites with V7-curves having Negative Differential Resistivity (NDR) 211]. Since JV-curves in semiconductors map to Vl-curves 
in superconductors, then NDC (for semiconductors) maps into NDR for superconductors. Here, N(S)-type shapes for semi- 
conductors correspond to S(N)-type for superconductors 3]. The Negative Differential Conductivity (NDC) found in ilO] is 
analogous to the Gunn effect in semiconductors, where electron-charge modulations lead to steps in j{E) in the NDC regime. 



-/Vint — N-u—Np, start to move and the velocity v becomes 
nonzero and grows with j (phase II in Fig. 2). With fur- 
ther increasing the current j, the driving force acting on 
a single vortex, Fd = j$o/c, overcomes the pinning force, 
and a significant fraction of the vortices start to move. 
This motion is uniform and very disordered. Phase III 
corresponds to such vortex-flow mode. At higher vortex 
velocities, the random motion of vortices becomes more 
ordered. Some vortices become pinned in commensurate 
rows while others move along vortex rows, which are in- 
commensurate with the underlying PAPS. Namely, when 
j exceeds a threshold value, only incommensurate vortex 
rows move and the vortex velocity exhibits a significant 
drop with the increase of the driving force (phase IV). 
Note that the phases III and IV have an analogy with 
the NDR behavior of electron motion in semiconductors 
with the increase of the voltage. Namely, increasing the 
applied force on the moving particles produces a lower 
velocity in them. At high current densities, the driving 
force completely overcomes the pinning force (phase V) 
and the curve v{j) tends to a Hnear one. 

The results obtained in Refs. [20 andl2ll prove that the 
V7-curves of superconductors with a square PAPS have 
a part with NDR of the S-type, since the electric field in 
the sample is related to the average vortex velocity by the 
well-known relation E = —v{j)B/c. Such type of NDR is 
usual for plasmas [y, 0| and semiconductors [l|, 0] giving 
rise to important instabilities of the uniform current flow 
known as pinch-effect and filamentary instability^ when a 
current flow breaks into filaments with lower and higher 



current densities [l|, [3, S 13, S fi 0| ■ ^^ superconductors 
we have only few examples of the S-type NDR in the 
samples with the specific weak links |10| . 

The described dynamical phases disappear in the case 
of very disordered pinning arrays and, consequently, the 
NDR of S-type in the VI -curves vanishes. Under realistic 
experimental conditions, the properties of the supercon- 
ductor in the flux flow regime are strongly affected by 
Joule heating since the current density j necessary to 
overcome the pinning force is high [5| . An increase of the 
sample temperature T due to Joule heat, jE, gives rise to 
a decrease of the pinning force, and the current density 
can drop down with the growth of the electric field. As a 
result, the V7-curve with an NDR of N-type (red dashed 
line in Fig. 2b) is commonly observed in superconductors 
for high current density [J, Q . The uniform state in sam- 
ples with NDR of N-type is also unstable [l|, 0] ; and a 
propagating resistive state boundary or the formation of 
resistive domains can destroy the uniform flux flow mode 
in the superconductor [1, d, Q . Increasing the pinning 
force and/or decreasing the thermal coupling of the sam- 
ple with its environment (which decreases the cooling 
rate of the sample), one can achieve a situation where 
the NDR of both N- and S-ty pe simultaneously coexist 
in the V7-curve (Fig. 2b) [21[. In this case, we predict 
remarkable flux flow instabilities. Note also that the ef- 
fect of thermal fluctuations on the flux flow regime in 
the superconductors with PAPS is somewhat analogous 
to the effect of positional disorder in the pinning sites. 
This effect should be taken into account, especially, at 



temperatures close to the critical temperature of the the 
superconductor T^. 

The effect of Joule heating on the VT-curve of super- 
conductors with PAPS was first outlined in Ref. 21. This 
gives rise to the coexistence of both types of instabili- 
ties, which is very unusual since non-linear devices are 
typically either N-type or N-type, but not both. Here 
we present a detailed analysis of this problem. In ad- 
dition, we study in detail the effect of disorder on this 
phenomenon. 

The paper is organized as follows. In Section II we for- 
mulate the model of the flux motion taking into account 
the effect of Joule heating. The heating manifests itself 
in thermal fluctuations of the vortices and a variation 
of the superconductor parameters due to the tempera- 
ture increase. Both of these effects are included in the 
model. In Section III the average velocity of the vor- 
tices, V, is found as a function of the applied current 
density j. This dependence is obtained by means of the 
molecular dynamics integration of the equations of mo- 
tion presented in Section II. As a result, we obtain the 
V7-curves of the sample with PAPS and study the ef- 
fects of temperature variation, thermal fluctuations, and 
positional disorder of the pinning sites on these curves. 
In Section IV an analytical criterion for the development 
of the filamentary instability in superconductors with an 
NDR of S-type is derived. In Section V we analyze the 
effect of the interplay between S-type and N-type NDR 
on the flux flow in superconductors with a PAPS. We 
argue that the coexistence of the NDR of both types can 
give rise to macroscopic non-uniform self-organized dy- 
namical structures in the flux flow regime. 




FIG. 1: The model: vortices driven by a Lorentz force pro- 
duced by an applied current in a superconductor with square 
array of pinning sites. The period of the pinning array is 
a — 2Ao. The pinning sites are shown by dark blue (dark 
grey) parabolic bars and by dark blue (dark grey) dots on the 
upper surface. The black dashed lines connecting the dots on 
the top surface are a guide to the eye. The vortices are shown 
by red-to-yellow (grey-to-light grey) tubes and by red (grey) 
larger spots on the top surface. The direction of the applied 
current j is indicated by a wide light blue (light grey) arrow, 
and the Lorentz force acting on the vortices Fa is shown by 
small red (grey) arrows. The direction of the external applied 
magnetic field B is shown by the brown (grey) arrow. 



unit length acting on the ith vortex due to the interac- 
tion with other vortices. The force per unit vortex length 
F^'' describes the interaction of the ith vortex with the 
pinning array. The term F^ arises due to the thermal 
fluctuation contribution to the force. As in standard ap- 
proaches Ff{t) is a random function of time t, obeying 
the correlation relations 



II. MODEL 

We describe the flux motion in a three-dimensional 
(3D) superconducting slab, infinite in the xy-plane, us- 
ing a 2D model (assuming no changes in the z-direction). 
This approach has also been used in the past, e.g., in 
Refs. ^, ^, ^, 25. We consider a sample with a square 
array of Np pinning sites interacting with N^ vortices re- 
lated to the magnetic field by B = Ny^Q. The magnetic 
field is perpendicular to the slab (see Fig. 1). The pe- 
riod of the regular array is a = 2Ao and we focus on the 
case when the magnetic field B is slightly higher than the 
first matching field i?^, that is N^ > Np. The vortices 
are driven by the Lorentz force, Fd = j$o/c, produced 
by the current flowing in the x direction (see Fig. 1). 
Thus, the horizontal axis of figures 2, 3, 4, and 5, refer 
to the driving force F^ or driving current j, since these 
arc proportional to each other. The overdampcd motion 
of the ith vortex is described by the equation 



Wi 



(1) 



where v^ is the velocity of ith vortex, rj — cr„_ffc2(r)$o/c^ 
is the fiux fiow viscosity, (T„ is the normal conductivity, 
and Hc2 is the upper critical field. F™ is the force per 



and 



{Fnt))t=0 (2) 



{F^{t)Fj{t'))t - 2TikBT5,,5{t - t'), (3) 



where ks is the Boltzmann constant, (...)( denotes a 
time average, 5ij is the Kronecker symbol, and 5{t) delta- 
function. 

We describe the vortex- vortex interaction by the usual 
expression for Abrikosov vortices 
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Ki 



A(r) 



(4) 



where A is the magnetic field penetration depth, Ki is 
the first order modified Bessel function, the summation 
is performed over the positions Tj of Ny vortices in the 
sample, and Vij — (r^ — rj)/\ri — rj\ is a unit vector in 
the direction of the force acting between the zth and jth 
vortices. 

The Np pinning sites (narrow indentations or "blind 
holes" which can accomodate a maximum of one vor- 
tex, for the vortex densities used in our calculations) are 

(p) 
located at positions r), . The pinning potentials are ap- 
proximated by parabolic wells. Then, the pinning force 



per unit length acting on the zth vortex can be written 
in the form 



Fp{T) 
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where r^ is the size of the elementary pinning poten- 
tial well, Fp(T) is the maximum pinning force, O is the 
Heaviside step function, Aq — X{T = 0), and r^-^' — 



{r^ 



ri^^)/|r. 



Ml 



is the unit vector in the direction 



of the elementary pinning force. In what follows, we es- 
timate the maximum pinning force as 



FpiT) 



Hl{T)e{T) 



(6) 



where ^(T) is the coherence length. 

The temperature dependence of the values entering 
our model is found using the Ginzburg-Landau approach. 
Therefore, Hc2{T) = ^o/2it^'^{T). The temperature de- 
pendence of the penetration depth A is approximated as 



A(r) = Ao ( 1 - ^ 



2\ -1/2 



(7) 



We assume that the Ginzburg-Landau ratio kql — A/f 
is independent of temperature. In this case, for the tem- 
perature dependence of the maximum pinning force, we 
have 



Fp(T) = Fp^[l-^ 



(8) 



where FpQ — Fpo(T = 0). We also assume that the nor- 
mal conductivity an is temperature independent. 

We simulate Eqs. H])-® using the molecular dynam- 
ics technique. Below we present the results for rectangu- 
lar cells with size 18 x 12 Aq. Periodic boundary condi- 
tions are imposed at the cell boundaries. First, we should 
prepare an initial state of our system. For this purpose 
we assume that the initial temperature of the system is 
high and that the vortex structure is in a liquid unpinned 
state. Then, we slowly decrease the temperature down 
to T = and vortices are captured by the pinning sites. 
When cooling down, vortices adjust themselves to min- 
imize their energy, simulating field-cooled experiments. 
Starting from this initial state, we increase the driving 
current and compute the average vortex velocity v{j), 
which is determined as 



^(j) = ^„~'E^'-*' 



(9) 



where x is the unit vector in the x direction. 

The equilibrium temperature distribution in the sam- 
ple can be found by solving the heat equation with Joule 
heating. The average power of this heating per unit vol- 
ume due to vortex motion is jE = jvB/c. The heat 



flux, q, removed from the sample boundaries by the ex- 
ternal coolant is described by the usual linear Kapitza 
law, q = Sho{T — To), where S is the sample surface, 
ho is the heat transfer coefficient, and To the ambient 
temperature. To simplify the procedure, we assume that 
the heat conductivity of the sample, k, is large, k:^ how, 
where w is the sample thickness. Under such a condition, 
the temperature in the sample is uniform and it can be 
found from the heat balance equation [5|: 



hoS{T-To) 



-JBV, 



(10) 



where V is the sample volume. Further, we shall assume 
that To < Tc and neglect To. 

We now introduce dimcnsionlcss variables. The dimcn- 
sionless current (which is equal to the dimensionless driv- 
ing force fd), dimensionless average vortex velocity T4, 
and dimensionless magnetic field induction &, are given 
by 



Jd= —, Kr = — , 0= — -. 
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The normalization values jo and vo are defined by 
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In dimensionless units the heat balance equation (jlOp re- 
lates the temperature and the dimensionless driving force 
as 

T 



— = Kth Vx fd b, 

J r. 



where 



i^th- 



jo Vq B^ V 

choTcS 



(13) 



is the ratio of the characteristic heat release to heat re- 
moval. In our calculations we used the values of the pa- 
rameters characteristic of high temperature superconduc- 
tors: Ao = 2000 A, KGL = 100, an = 10^^ s-\ V/S = 
1000 A, Tc = 90 K, B^ = 500 G, and ho = 1 W/cm^K. 
In this case, we find that Kth = 0.05 — 0.06 and Fpo is 
of the order of Tq ~ <I>oio/c. In the simulations, we used 
Kth = 0.0525 and Fpo = 2To . 



III. SIMULATION RESULTS 

A. Effect of heating: V7-curve of N- and S-type 

The calculated dependence of v versus the current den- 
sity j in the absence of the Joule heating is shown in 
Fig. 2a, in dimensionless units. The shape of this curve 
is similar to that as found in Ref . [20J and [2l|: there ex- 
ist five dynamical vortex phases described in the intro- 
duction and a pronounced hysteresis for increasing and 
decreasing current regimes. 



A significant effect of tfie heating is observed if the cur- 
rent density exceeds some threshold value, j > 3.75 for 
the case shown in Fig. 2b. In particular, the jump from 
phase II to phase III becomes larger, and an abrupt tran- 
sition occurs between regimes IV and V, compared to the 
non-heating case shown in Fig. 2a. The most important 
feature related to Joule heating in the high current range 
is the appearance of hysteresis in some regions of phases 
IV and V. For decreasing current, the overheated vortex 
lattice keeps moving as a whole at lower currents than 
the "cold" one (for increasing j). The transition part of 
the v{j) curve from phase IV to phase V is shown by 
the dashed line in Fig. 2b. This part of the V7-curve is 
unstable for a given drive force and could be only found 
for a fixed voltage Q . 

As a result, we obtain a new complex NDR of a hybrid 
nature with both N- and 5-type instabilities, which is 
very unusual for any media, especially for superconduc- 
tors. Each type of NDR is characterized by its specific 
instabilities |l|, [3, i, i, i, i, 0]. Thus, the obtained VI- 
curve is characterized by two kinds of instabilities. For 
example, if the current density exceeds the value j ~ 3.5 
(point A in Fig. 2b), the uniform current flow becomes 
unstable and a filamentary instability [l|, |j] occurs. Due 
to this instability, the current flow breaks into filaments 
with different supercurrent density, some with lower cur- 
rent JB (state B) and others with higher current jc (state 
C). However, the state C is, in its turn, unstable and de- 
cays. The corresponding stable states are on the lower 
(E) and on the upper (D) V7-curve branches. That is, 
the filament breaks into domains with higher and lower 
value of vortex flow speed v; in other words, with higher 
and lower electric field. The stability and evolution of 
such a complicated structure is an open question (see 
also Section V). 

In Fig. 3, the average vortex velocity v{j) is shown 
for different radii of the pinning sites. The shape of the 
function v{j) only slightly changes for the radii in the 
range r^ = 0.20Ao to 0.25Ao (Figs. 3a, b). However, 
for radii smaller than a certain value, phase IV in the 
reverse branch (i.e., when decreasing the driving current 
j) disappears (Fig. 3c) since the overheated vortex lattice 
cannot adjust itself to the pinning array and turns to the 
disordered motion in phase HI. 



B. Effect of disorder 
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FIG. 2: The average vortex velocity v (x E versus current j 
for B/B^ = 1.074, rp = 0.21Ao and Fpo/Fo = 2 for increas- 
ing (shown by red (black) solid circles) and for decreasing 
(blue open squares) j (see, Ref. [2l|). (a) No heating effect 
is taken into account. The regions corresponding to differ- 
ent phases are indicated by the roman numerals from I to V 
(as in Ref. [20l. [21[|). (b) The sample heats up due to vortex 
motion. For small values of the drive j (up to j ~ 3) the 
effect of heating is negligible. For j ~ 3.2, a jump from phase 
II to phase III occurs. As a result of the heating, a transi- 
tion occurs abruptly at j ~ 4.3, from regime IV to regime V, 
where the vortex lattice is entirely unpinned and moves as a 
whole. A hysteresis now appears in region V: when decreas- 
ing the current j down to the value at which the jump from 
phase IV to phase V occurred, when the driving increased, 
the overheated vortex lattice keeps moving as a whole. As 
a result, we obtain a complicated A^- and 5- type Vl-curve. 
State A is unstable and the sample divides into filaments in 
states B and C. State C is also unstable. The correspond- 
ing stable states are on the lower (point E) and on the upper 
(point D) V7-curve branches. Point j ~ 4.5 corresponds to 
the normal transition and at T > Tc we have the usual Ohmic 
conductivity. 



Let us study the effects of disorder on the NDR. Small 
disorder can be effectively introduced to the system by 
increasing the radius, Vp, of the pinning sites. So vortices 
acquire an additional degree of freedom and can move 
inside the pinning sites. The function v{j) is shown in 
Fig. 4 for larger pinning site radii, Vp = 0.35Ao (Fig. 4a), 
and rp — 0.45Ao (Fig. 4b). In case of larger radii, phase 
HI disappears, and the motion of interstitial vortices 
(phase II) transforms directly to the ID incommensu- 
rate vortex motion (phase IV) (Fig. 4a). However, the 



robust hysteresis related to heating remains. For larger 
radii of the pinning sites, phase II, related to the motion 
of interstitial vortices, disappears since all the vortices 
are pinned for weak enough drives (Fig. 4b). 

To model disorder related to a distortion of the reg- 
ular PAPS, we introduced small random displacements 
for each pinning site. Specifically, for the displacement 
of each pinning site we consider a random angle ctran 
(0 < aran < Stt) and a random radius Tian (0 < rj-an < 
r^fj'^, measured in units of a/2, where a is a period of the 
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FIG. 3: The average vortex velocity v{j) for B/B^ = 1.074, 
Fpo/Fo = 2.0 for increasing (red/black solid circles) and 
for decreasing (blue open squares) j and different radii of 
the pinning sites: (a) rp = 0.23Ao; (b) rp = 0.20Ao; (c) 
Vp = O.IQAq. The function v{j) only slightly changes for radii 
from rp — 0.20Ao to 0.25Ao. For smaller Vp, phase IV in the 
reverse branch disappears (c) since the overheated vortex lat- 
tice cannot adjust itself to the pinning array. 



(regular) pinning array) for each pinning displacement. 

The corresponding values of v{j) are shown in Fig. 5 
for different amounts of disorder. Note that small disor- 
der does not appreciably influence v{j) (Fig. 5a). How- 
ever, for r^^ > 0.05(a/2) (Fig. 5b), phase IV disap- 
pears in the reverse branch of v{j). For r™^'^ > 0.1(a/2) 
(Fig. 5c), phase IV is lost in both branches; only some 
reminiscent features remain, which disappear for larger 
d^T (Figs- 5d, e). Finally, at full disorder r™^'' = a/2 
(Fig. 5f), only phases I, III and V remain. However, even 
in this case there is a weak hysteresis related to heating 
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FIG. 4: The average vortex velocity v{j) for increasing (red 
(black) solid circles) and for decreasing (blue open squares) j 
for large pinning sites radii: (a) rp = 0.35Ao; (b) rp — 0.45Ao. 
Other parameters are the same as in Fig. 2. An increase of 
rp increases the disorder in the system because vortices can 
then move inside the pinning sites. Phase III disappears in 
(a). However, the hysteresis related to the heating remains. 
For larger rp (b), phase II disappears since all the vortices are 
pinned. 



(inset to Fig. 5f), observed in experiments with random 
pinning Q. 



IV. FILAMENTARY INSTABILITY 

It is well-known that the uniform current and electric 
field distributions are unstable under definite conditions 
if the sample V7-curve has parts with NDR [l|, 0, S 0, 
H) S Q- Let us assume that the sample with the VI- 
curve shown in Fig. 2b is in a current-biased regime. If 
the driving current exceeds the value j w 3.75, the uni- 
form current flow with the current density j becomes 
unstable with respect to the so-called filamentary insta- 
bility [l|, 0. Thus, the current flow in the sample breaks 
up into stripes or filaments with two different alternating 
current densities. This process is illustrated in Fig. lb: 
the sample in state A with current density ja breaks into 
filaments with lower jb (state B) and higher jc (state C) 
current densities. 

To study the process in more detail, we consider a sam- 
ple connected to a standard electrical circuit, Fig. 6. The 
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FIG. 5: The average vortex velocity u as a function of 
the driving current j for B/B^ = 1.074, r^ = 0.2Ao and 
Fpo/Fo — 2.0 for increasing (red (black) solid circles) and for 
decreasing (blue open triangles) j for different amounts of dis- 
order (displacement of the centers of the pinning sites from 
their regular positions) in the system: (a) dran ~ 0.01(a/2); 
(b) dran = 0.05(a/2); (c) d,an = 0.1(a/2); (d) d^an = 0.2(a/2); 
(e) dran = 0.5(a/2); (f) dran ~ ffl/2. The function v{j) does 
not appreciably change for small amount of disorder (a) . For 
dran > 0.05(a/2) (b), phase IV disappears in the reverse 
branch. For dran > 0.1 (a/2) (c), phase IV is lost in both 
branches; only some reminiscent features remain, which dis- 
appear for larger dran (d, e). At maximal disorder dran ~ a,/2 
(f), only phases I, III and V remain. However, there is a weak 
hysteresis related to heating (inset to (f)). 
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FIG. 6: The electrical circuit. The current filaments are 
shown schematically by light grey lines [23|; I is the sample 
length, w is the sample half-width, Au; is the filament width, 
the boundary conditions are stated at the sample edges y = 



circuit equation is 



dl 



RI + IE = U, 



(14) 



where / is the current in the circuit, L and R are the 
circuit inductance and resistance, I is the sample length, 
E is the electric field in the sample, U is the voltage at 
the circuit terminals, which is assumed to be constant, 
and j = I /A, where A is the sample cross-section. Us- 
ing Eq. p^ and Maxwell equations we can write the 
equations describing the development of the small per- 
turbations of electromagnetic field 6E, JB, and current 
density Sj in the form 

LA-^^{Sj) + RA{S]) + lSE = 0, 

Vx.E.-ii^, (15, 

c ot 

VxSB = —{dj). 
c 

Here 

all quantities are assumed averaged over a volume in- 
cluding a large number of vortices, Sj and 6E denote the 
average values over the sample cross-section. Below we 
assume for simplicity that the sample has zero demagne- 
tization factor. 

We shall seek the solution to Eqs. ((T5| in the standard 
form: 5E, SB ex exp(Ai/io); where A is the value to be 
found, and fo — L/ R is the circuit relaxation time. 

An instability develops if Re (A) > 0. In general we 
should add to Eqs. P3)) the equation for the small tem- 
perature perturbations but here we study the filamentary 
instability for which the temperature rise is not of crucial 



importance. We also neglect the self-field effect and as- 
sume that the background magnetic field in the sample is 
uniform. To find the filamentary instability criterion we 
can consider the perturbations depending on the y coor- 
dinate only [l|, y] ■ In such a geometry, the perturbation 
of the electric field has only the x component, while (5B 
has only the z component. Using 6B = ct^SE' /Xw, we 
find from Eqs. P^ 



{\ + \)S']+p-^5E = 0, 



SE" - PSE' 



to 



SE = 0, 



(16) 
(17) 



where prime means differentiation over the dimensionless 
coordinate y/w, w is the sample half-width, pc — RA/l, 
and 



ts 



c2 ae 



/3 



A'Kw dj 



are the sample VT-curve parameters, which are positive 
or negative depending on the relevant part of the VI- 
curve at given background fields E and B. Note that 
the value of \ts\ is the characteristic time of the magnetic 
field relaxation in the sample. 

The differential equation (I17p requires two bound- 
ary conditions. We obtain the first one assuming that 
the applied magnetic field B is constant. This means 
that 5B{1) = -SB{-1) or 5E'{1) = -SE'{-1). From 
the MaxweU equations ^Bj we get 5B{1) - SB{-1) = 
8ttwSj/c. Using these 6B^s, and substituting dj from 
Eq. (|16p . we obtain the second boundary condition in 
the form 



<5^'(1) 



jXSE 
(A+T)' 



where 



7 



Attw 



c?AL' 
The solution of Eq. ((T71) reads 

5E ^ Ci exp(piy/w) -t- C2 exp(p2y/fi'), 

where Ci are constants and 



(18) 



P /?2 Ai, 

Substituting Eq. (fTSjl to the boundary conditions we ob- 
tain a set of uniform linear equations for the constants 
Ci^2- The non-trivial solution of this equation set exists 
if its determinant is zero. Thus we find the equation for 
the eigenvalue spectrum A in the form 



Pi 



P2 



P2 



Pi 



7A 1 


(A+1)P2J 


7A 1 


(A + l)piJ 



coshpi sinhp2 = 



sinhpi coshp2- 



(19) 



In the simplest case of small \dj/dB\, when 



-P2 = \/Xts/tQ, 



Pi 



the solution of Eq. (fTO|) can be readily found explicitly 
with an accuracy up to \0^\ 



A = 



ts d] 



It follows from the last relation that the instability occurs 
only at the V7-curve branch with NDR when 

ts ex -^ < 
dE 

and, moreover, the drop of the voltage should be large 
enough, 



dE 



dj 



> pc 



(21) 



In this case pi.2 are purely imaginary numbers and the 
solution of Eqs. (14) is periodic in the j/-direction. The 
characteristic spatial period of the arising current fila- 
ment structure is of the order of 



Ai 



biWI VW+tJhl' 



(22) 



In a very unstable state, \jts/to\ ^ 1, the filament width 
is small. Aw <^ w. The characteristic instability build- 
up time becomes io/A. Thus, a sample with an S-type 
NDR in VT-curve divides itself into small filaments with 
different current densities (in different dynamic flux flow 
phases III and IV) if the resistance and inductance of the 
external circuit are restricted by inequalities 



R«^ 



dE 



dj 



L <c AttI/c^A . 



(23) 



The obtained results are valid if 



_^dE_ 

'" a7 



> 



Airw dj 
c dB 



According to Eq. ([2T|) . the left hand side of the last 
inequality should be higher than unity, while the right 
hand side is much smaller than 1 for the parameter range 
studied in the previous sections if the sample half-width 
w < 1 mm. Note, however, that taking into account the 
magnetic field dependence of the VT-curve gives rise to 
some increase in its stability. 



V. INTERPLAY BETWEEN N-TYPE AND 
S-TYPE INSTABILITIES 

In the stationary inhomogeneous state that arises after 
the development of the filamentary instability, the elec- 
tric field should be uniform over the sample. The part, p, 
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FIG. 7: Schematic V7-curves (red (black) solid line and green 
(dark grey) dashed line) of the superconductor for two differ- 
ent values of the hysteresis due to overheating [2j| . The more 
pronounced hysteresis loop (shown by the green (dark grey) 
dashed line) corresponds to a larger value of the characteristic 
heat release. 



of the filament with the higher current jc and the part, 
1 — p, with the lower current jb (see Fig. 2b) are defined 
by an evident condition ja = pjc + (1 ^ p)JB, that is, 



P 



Uc-Ja) 
Uc - 3b)' 



For the case shown in Fig. lb we find an estimate p ~ 0.6. 

As shown in Fig. 4, sufficiently high disorder destroys 
the NDR of S-type in the F/-curve, but the NDR of N- 
type can still exist in this case due to sample heating. 
Such an unstable regime has been thoroughly studied for 
superconductors 0,0,0, and we do not discuss it here. 

A richer dynamics can be observed if the VT-curve has 
NDR parts of both N- and S-types, Figs. 2b and 3. In 
this case the filaments with higher current density jc 
are unstable if the system is far from the voltage-biased 
regime [l|, Q]. The instability of the filament with N- 
type V7-curve should switch the filaments into state D 
(Fig. 2b) with high resistivity or to the formation of the 
domain structure with higher, D, and lower, E, resistiv- 
ity [J| . However, any possible decay of the unstable state 
C gives rise to a non-uniform electric field distribution in 
the sample and, as a result, to non-zero dB/dt. Thus, 
the state that appears after the instability develops is not 
stationary but rather a dynamic one. 

To clarify the situation, we consider two possible VI- 
curves shown in Fig. 7 (one red and another, labelled 
by " 2" , with a green branch) and assume that the to- 
tal current value in the circuit is fixed, / = Aja- In a 
more general case, after the decay of the unstable state 
C, the filaments with higher current density will be over- 
heated and transit to the higher resistive state D. In this 
state D the flux lines move fast, which, along with the 
temperature increase due to the thermal conductivity, 
gives rise to the acceleration of the fiux flow in the lower- 
current filaments. As a result, the high-resistivity over- 



heated state moves from point D to a lower electric field 
range, while the low-resistivity state (point B) moves to 
a higher electric field range (point A). If we assume that 
the system has a VT-curve of the type 2 (green (dark 
grey) dashed curve in Fig. 7), then the high-current fila- 
ments in D move to the state A' with current density ja 
and a higher electric field than in the initial state A. The 
filaments in point B move to point A and jump to point 
A'. As a result, a new stable uniform state with j — ja 
appears. However, if the VT-curve has the form shown 
by the red curve 1, the stable stationary point with the 
current density ja does not exist. In this case, the high 
electric field state moves from point D to the marginal 
stability point F and then falls down to a lower branch of 
the VT-curve curve (point F'). In this state F' the elec- 
tric field is lower than in the state B with lower current. 
As result, the state moves from F to F' and then to point 
A, which is the only uniform state corresponding to the 
fixed current value ja- However, this state is unstable, 
and the cycle 



C^D 



F'^A 



is repeated (branch 1 in red in Fig. 7). Another branch 
(branch 2 with the green segment) of this cycle involves 
the stable filament state B and is as follows: 



A^B 



A. 



Such a dynamic state has a non-stationary pattern of re- 
sistive domains coexisting and intertwinned with current 
filaments. The specific form of these patterns, and their 
dynamics can be very complex, and requires investiga- 
tions beyond the scope of this study. 



VI. CONCLUSIONS 

The influence of temperature on the dynamic phases 
and current-voltage characteristic of superconductors 
with periodic pinning array was investigated here. It is 
demonstrated that this effect can change the VT-curve 
drastically. For a range of values of the pinning array 
parameters and heat transfer characteristics it is possi- 
ble to obtain the VT-curves with a negative differential 
resistivity of (i) either N- or S-type, (ii) VT-curves with 
both types of NDR, or VT-curves without any NDR parts. 
The uniform flux flow is unstable if the VT-curve has a 
part with NDR. The formation of resistive domain struc- 
tures and/or propagating the resistive state through the 
sample is a characteristic of VT-curves of N-type, while 
a filamentary instability with sample regions (filaments) 
having different current densities is a characteristic of 
VT-curves with NDR of S-type. Much more complex 
regimes can be expected in the case of VT-curves with 
NDR parts of both types. In this case the possibility of 
arising dynamical non-uniform regimes is argued. 
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